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非 线 性 微分 动力 系统 稳定 域 计算 的 波形 松弛 方法 * 


菌 小 林 
(西安 交通 大 学 理学 院 ， 西 安 710049) 


d E 非 线性 微分 动力 系统 稳定 域 计 算是 在 许多 领域 具有 实际 应 用 的 问题 。 本 文 对 非 线 性 微分 动力 系统 
稳定 域 的 计算 方法 进行 了 总 结 ， 通 过 对 稳定 域 边界 流 形 的 全 面 分 析 ， 提 出 了 用 波形 松弛 方法 计算 
稳定 域 边界 流 形 的 思想 ， 给 出 了 计算 稳定 域 边界 流 形 的 波形 松弛 算法 。 第 一 步 ， 找 出 微分 动力 系 
统 的 所 有 平衡 点 ， 确 定 渐 近 稳 定 平 衡 点 稳定 域 边界 上 的 平衡 点 ， 第 二 步 ， 用 微分 动力 系统 的 反方 
向 系统 确定 原 系统 渐 近 稳定 平衡 点 稳定 域 边界 上 平衡 点 的 稳定 流 形 ， 第 三 步 ， 渐 近 稳 定 平 衡 点 稳 
定 域 的 边界 是 由 边界 上 平衡 点 和 该 平衡 点 稳定 流 形 的 并 集 构成 。 最 后 用 例子 来 说 明 波形 松弛 方法 
在 计算 稳定 域 边界 流 形 的 有 效 性 。 
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非 线性 微分 动力 系统 平衡 点 的 稳定 域 ( 或 称 为 吸引 域 ) 在 诸如 电力 系统 ， 经 济 系 统 ， 生 
物 系 统 等 方面 都 有 非常 重要 的 应 用 。 从 应 用 的 观点 来 看 希望 得 到 最 大 的 稳定 域 ， 这 样 在 
电力 系统 的 设计 叫 及 相关 参数 的 选取 中 可 以 节省 由 于 不 必要 的 过 保护 设计 而 产生 的 浪费 。 
同样 在 经 济 系统 、 化 学 反应 器 、 生 物 系统 四 等 方面 可 以 节省 大 量 的 人 力 物 力 ， 实 现 最 大 
的 经 济 效益 。 确 定 稳 定 域 的 方法 有 Lyapunov 函数 方法 、 拓 扑 动力 系统 方法 和 其 它 数值 方 
法 等 ， 此 方面 比较 早 的 综述 文章 有 Genesio，Tartaglia，Vicinolal 和 Varaiya, Wu, Chen, 
他 们 总 结 了 二 十 世纪 八 十 年 代 以 前 稳定 域 的 研究 工作 ， 给 出 了 计算 稳定 域 的 方法 和 所 
取得 的 成 果 ， 同 时 基于 轨 线 的 拓扑 结构 提出 了 “ 轨 线 逆转 ”方法 。Lyapunov 函数 方法 就 
是 通过 对 所 研究 的 微分 系统 构造 相应 的 Lyapunov 函数 ， 应 用 稳定 性 方面 的 结论 得 到 系统 
渐 近 稳定 平衡 点 的 稳定 域 ， 这 种 稳定 域 具 有 较 大 的 保守 性 ， 通 常 是 稳定 域 的 子 域 ， 如 
通过 对 Lyapunov 函数 进行 逐次 迭代 来 得 到 稳定 域 的 一 部 分 以 及 通过 构造 最 大 Lyapunov PS 
数 方法 回来 得 到 稳定 域 等 。 最 近 ，Chesil6? 提出 了 用 一 系列 连续 Lyapunov 函 数 而 不 是 一 
个 Lyapunov 函数 来 估计 多 项 式微 分 系统 的 稳定 域 。 拓 扑 动力 系统 方法 主要 是 基于 微分 拓 
扑 和 动力 系统 的 基本 理论 ， 在 相 空间 对 轨 线 的 全 局 结构 进行 分 析 ， 进 而 给 出 稳定 域 的 计算 
方法 ，Chiang, Hirsch 和 Wulsl，Chiang #0 Thorpl?!, Lee Chiang 以 及 Zaborszky，Huang， 
Zheng fll Leung!!! 等 从 拓扑 动力 学 的 观点 出 发 对 稳定 域 进行 了 深入 研究 ， 取 得 了 显著 的 结果 。 
计算 稳定 域 还 有 其 他 方法 ， 如 能 量 函 数 方法 、 数 值 方法 等 。 

波形 松弛 方法 是 微分 方程 数值 解法 ， 它 将 整个 微分 系统 分 成 若 于 个 子 系统 ， 各 自 独 立 进行 
仿真 计算 ， 子 系统 之 间 的 耦合 关系 通过 子 系统 中 的 一 些 波形 信号 的 相互 传递 来 实现 ， 在 系统 之 
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间 联 系 相 对 较 弱 的 情况 下 ， 比 较 容 易 得 到 收敛 结果 。 

本 文 对 微分 动力 系统 的 稳定 域 进行 了 综合 分 析 ， 得 到 微分 动力 系统 一 个 渐 近 稳定 平衡 点 稳 
定 域 的 边界 是 由 闭 轨 线 或 由 不 稳定 平衡 点 以 及 该 不 稳定 平衡 点 的 稳定 或 不 稳定 流 形 构成 。 基 于 
这 些 结果 给 出 了 寻找 稳定 域 边 界 流 形 的 波形 松弛 数值 算法 ， 并 用 例子 说 明 波形 松弛 方法 在 计算 
稳定 域 边 界 流 形 的 有 效 性 。 


2 动力 系统 基本 概念 与 波形 松 驰 算法 


对 于 非 线 性 自治 微分 动力 系统 
dz 
T fo. o) 
dz 


其 中 zeMCR" f : M 一 M， 我 们 假设 向 量 场 是 连续 函数 并 且 对 其 变量 具有 一 阶 连续 偏 
导数 ， 以 保证 系统 (1) 和 (2) 解 的 存在 唯一 性 ， 在 相同 时 间 变 化 情况 下 ， 系 统 (1) 和 (2) 的 轨 线 
走向 方向 相反 ， 因 此 称 系统 (2) 为 系统 (1) 的 反方 向 系统 。 

方程 


f(z) =0 (3) 


的 一 个 解 称 为 系统 (1) 的 一 个 平衡 点 ， 我 们 用 集合 已 = {x : f(x) = 0} 来 表示 系统 (1) 的 平衡 点 
集合 ， 巨 中 既 有 孤立 点 也 可 能 含有 连续 可 微 闭 轨 线 。 

在 局 部 坐标 系 下 ， 如 果 在 平衡 点 ze 处 的 雅 可 比 矩 阵 .fF 的 特征 值 没 有 零 实 部 ， 则 称 平 
衡 点 ze 为 系统 (1) 的 双 曲 平衡 点 。 对 一 个 双 曲 平衡 点 ze， 可 以 把 系统 (1) 在 点 ze 处 的 切 空 
间 7Tz.(M) 唯 一 地 分 解 成 两 个 子 空间 和 E* 的 直 和 ， 即 Bs @ E*， 每 个 子 空间 在 线性 算 
子 Jz。f 下 是 不 变 的。 对 于 ESK, Je f 的 特征 值 具有 负 实 部 ， 对 于 ERKA, Ju, f 的 特征 值 
RAER. WEWEKA ns WE RERA no TZR Es 和 Er 可 表示 如 下 

稳定 的 子 空间 Es = span{ful uv2 ... ,v"*), 
不 稳定 的 子 空间 E"-span(wl,w?,... w}. 
这 里 ,v2,…. ,vm 是 ms 个 ( 正 交 化 ) 特征 向 量 ， 它 们 的 特征 值 具 有 负 实 部 ， 而 让 1 w, o, w 
E nu, A (EZE) 特征 向 量 ， 它 们 的 特征 值 具 有 正 实 部 ， 显 然 mn。 + mu = n. 
对 系统 (1)， 称 平衡 点 ze 为 wu 型 。 一 个 零 型 的 平衡 点 称 为 汇 ， 一 个 ?型 的 平衡 点 称 为 源 ， 


其 它 型 的 平衡 点 称 为 鞍点 ， 显 然 汇 是 渐 近 稳定 的 平衡 点 ， 源 和 鞍点 是 不 稳定 的 平衡 点 。 
设 是 一 个 双 曲 平衡 点 ， 它 的 稳定 流 形 Ws() 和 不 稳定 流 形 W*(2) 定义 如 下 


W*(3)— {z € M : (z) 5$, 4 t coo], (4) 

W*(£)- (re M :&(z) ê, 4 t —oo). (5) 
类 似 地 ， 一 个 双 曲 型 闭 轨 的 稳定 流 形 Ws(y) 和 不 稳定 流 形 W” (y) 定义 如 下 

W*(y) 2(z€ M:&i(z) ^5 y, 3 t coo], (6) 


W*(y)—-íreM:&,((z)-«, 当 t5 —oo]. (7) 
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H F Ai D) 临界 元 素 的 稳定 流 形 对 应 于 理 _:(:) 临界 元 素 的 不 稳定 流 形 ， 通 过 这 种 性 质 可 
以 把 稳定 流 形 的 性 质 转换 为 不 稳定 流 形 的 性 质 。 对 于 一 个 集 SCR"， 如 果 从 SS 中 出 发 的 
微分 动力 系统 (1) 的 每 一 条 轨 线 对 任何 时 刻 t 始 终 属于 5， 则 称 5 为 (1) 的 不 变 集 。 显 然 集 
tr W*()50W'(-) 8825 (1) 的 不 变 集 。 

轨 线 的 无 穷 远 性 质 可 以 用 它 的 w- 极 限 集 和 a- 极 限 集 来 研究 。 如 果 存 在 展 中 的 时 间 序 
Itip 3t; 一 +oo 时 ， 使 得 


y= ,m, Bt (z), 


JUR y 属于 z 的 w- 极 限 集 ( 记 为 w(z))。 类 似 地 ，a- 极 限 集 a(z) Hti 一 -co 来 定义 。 可 以 证 明 
这 些 极限 集 是 集合 M 的 闭 不 变 子 集 。 例 如 ， 一 个 平衡 点 是 它 自己 的 w- 极 限 集 ， 轨 线 的 w- 极 
限 集 在 它 的 稳定 流 形 内 ， 轨 线 的 a- 极 限 集 属于 它 的 不 稳定 流 形 。 一 个 闭 轨 Y 既 是 Y 上 任何 点 
的 w- 极 限 集 ， 也 是 7 上 任何 点 的 a- 极限 集 。 

横 截 性 的 概念 是 动力 系统 的 基本 概念 ， 如 果 4, BJÉ M 中 的 单 射 浸入 流 形 ， 且 满足 在 每 一 
点 TE ANB 处 ，A 和 B 的 切 空间 张 成 z 处 MM 的 切 空间 ， 即 


T,(A) +T,(B) = TM), ze AnB, 


RA, BB 完全 不 相交 ， 则 称 4 各 满足 模 截 性 条 件 。 

一 个 双 曲 平衡 点 的 重要 特性 之 一 就 是 它 的 稳定 和 不 稳定 流 形 在 $ 处 模 截 相交 。 横 截 相交 
的 重要 性 是 由 于 在 向 量 场 的 扰动 下 它 保 持 不 变 。 

对 系统 (1) 和 系统 (2) 作 如 下 假设 : 

AD 稳定 边界 上 的 所 有 平衡 点 是 双 曲 型; 

A2) 在 稳定 边界 上 所 有 平衡 点 的 稳定 和 不 稳定 流 形 满足 横 截 相交 条 件 ; 

A3) 稳定 边界 上 每 条 轨 线 当 t 一 co 时 趋 于 平衡 点 之 一 。 

考虑 微分 方程 的 初 值 问题 

[e ft,z), 


do (8) 
z(to) = To, te [to, T], 


HP f: [to, T] x R” >R”, zx € R”, RB FRERE HER. 
微分 动力 系统 (8) 一 般 的 波形 松弛 迭代 格式 为 
ic o0 2 F(t, g tV (t), PAG] (£)), 
g 1 (to) —2z9, t€ [to, T], 


JC Ab ER CP AE 
F :lto,T) x Rx R” R^, F(t,z,z)-— f(t, x). 


HOSCE RS 3877 OR VE AARAA ACA VE E BEER ATE. Jiangl 2:13! fl Gander, 
Ruehli!4 等 用 波形 松弛 方法 研究 了 微分 动力 系统 初 值 问 题 ， 给 出 了 和 迭代 格式 收敛 的 充分 性 条 
件 ， 用 波形 松弛 方法 、 单 调 波形 松弛 方法 和 Krylov 子 空间 波形 松弛 方法 等 对 微分 动力 系统 初 
值 问题 以 及 周期 边 值 问题 进行 了 研究 ， 得 到 了 选 代 序列 收敛 及 单调 收敛 的 充分 性 条 件 ， 这 些 充 
分 性 条 件 保 证 本 文 用 波形 松弛 方法 计算 微分 动力 系统 一 个 渐 近 稳定 平衡 点 稳定 域 边 界 流 形 算法 
的 收敛 性 。 
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3 ”确定 稳定 域 边 界 流 形 的 波形 松弛 算法 


根据 渐 近 稳定 平衡 点 稳定 域 边界 流 形 的 结构 特征 ， 在 假设 条 件 Al) 到 A3) 成立 情况 下 ， 有 
如 下 计算 稳定 域 边界 的 波形 松弛 算法 。 

算法 ”确定 渐 近 稳定 平衡 点 zs 的 稳定 域 边界 94(z。)。 

第 1 步 : 找 出 系统 (1) 的 所 有 平衡 点 ， 确 定 渐 近 稳 定 平 衡 点 x, 稳定 域 边界 上 的 平衡 点 ; 

第 2 步 : 用 系统 (2) 确定 系统 (1) 渐 近 稳定 平衡 点 z。 稳定 域 边界 上 这 些 平衡 点 的 稳定 流 形 ; 

第 3 步 : 渐 近 稳定 平衡 点 zs 稳定 域 的 边界 是 由 第 2 步 确 定 的 平衡 点 和 该 平衡 点 稳定 流 形 的 
并 集 构 成 。 

算法 中 第 1 步 包 括 找 出 f(z) = 0 的 所 有 解 。 第 2 步 可 以 用 波形 松弛 方法 来 完成 ， 过 程 如 
T: 

首先 ， 确 定 系 统 (1) 渐 近 稳定 平衡 点 zs 稳定 域 边界 上 的 平衡 点 。 
i) 找 出 系统 (1) 的 所 有 平衡 点 , 零 型 平衡 点 记 为 zs， 否 则 记 为 2;， 并 求 出 系统 (1) 在 平衡 
点 你 处 的 雅 可 比 矩 阵 ; 

i) 找 出 雅 可 比 矩 阵 具有 正 实 部 特征 值 且 具 有 单位 长 度 的 不 稳定 特征 向 量 ， 记 为 y; 

ii) 找 出 这 些 具有 单位 长 度 的 不 稳定 特征 向 量 y; 和 平衡 点 3; 的 e- 球 边界 的 每 一 个 交 截 ( 交 
BRAE ĉi 十 eyi FO ĉi 一 Eyi): 

iv) 从 每 个 交 截 点 开始 对 系统 (2) 用 波形 松弛 方法 计算 积分 轨 线 到 某 一 步 ， 如 果 轨 线 
保留 在 这 个 e- 球 内 ， 则 进行 下 一 步 ， 否则 ， 用 Qe 代替 e， 相 应 的 交 截 点 为 2; aey Ie 
处 0 < a < 1， 重复 这 一 步 ; 

v) ”从 这 些 交 截 点 开始 用 波形 松弛 方法 计算 积分 轨 线 ; 

vi) 重复 十) 到 v) 步 。 如 果 任 何 轨 线 趋 于 zs， 则 平衡 点 在 稳定 域 边界 上 。 

RR, MERR (1) 渐 近 稳定 平衡 点 z。 稳定 域 边界 上 平衡 点 的 稳定 流 形 。 

对 平面 系统 ， 稳 定 域 边界 上 平衡 点 的 类 型 是 一 个 鞍点 或 两 个 源 。 在 这 种 情 型 下 ，1- 型 平衡 
点 的 稳定 流 形 是 一 维 ， 它 可 以 用 波形 松弛 方法 确定 如 下 : 

a) ” 找 出 在 平衡 点 多 处 雅 可 比 和 矩阵 单位 长 度 的 稳定 特征 向 量 y; 

b) 找 出 这 个 单位 长 度 的 稳定 特征 向 量 y 与 平衡 点 人 的 e- 球 边界 的 交 截 ( 交 截 点 为 全 十 
ey FI ê — ey); 

c) 从 这 些 交 截 点 开始 对 系统 (2) 用 波形 松弛 方法 计算 向 量 场 的 积分 轨 线 到 某 一 步 。 如 果 轨 
线 保留 在 这 个 e- 球 内 ， 则 进行 下 一 步 ， 相 反 ， 我 们 用 ae 代替 =， 相应 的 交 截 点 为 全 十 aey， 此 
处 0 < a < 1， 重 复 这 一 步 ; 

d) ”从 这 些 交 截 点 开始 对 系统 (2) 用 波形 松弛 方法 计算 积分 轨 线 ; 

e) ”所 得 轨 线 是 系统 (2) 平衡 点 的 不 稳定 流 形 。 

对 高 维系 统 ， 计 算 过程 类 似 于 上 述 过 程 ， 计 算 结 果 能 提供 渐 近 稳定 平衡 点 稳定 域 边界 的 稳 
定 流 形 轨 线 集 。 

注 ”如 果 系 统 (1) 只 有 一 个 渐 近 稳定 平衡 点 ， 则 在 该 平衡 点 较 远 处 取 几 个 点 ， 从 这 些 点 处 开 
始 用 波形 松弛 方法 计算 积分 轨 线 ， 可 以 得 到 渐 近 稳定 平衡 点 稳定 域 边界 。 


4 数值 算 例 
例 1 考虑 如 下 系统 四 
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d 
E = 一 271 + Z172， 
(10) 
dzz = —T2 + T1T 
dt 2 1-22; 
dz, 
do 2zi 一 zir2, 
(11) 
dia = T2 — TIT 
dt = 142; 


容易 得 到 系统 (10) 的 两 个 平衡 点 (0,0) 和 (1,2)， 其 中 (0,0) 是 渐 近 稳定 平衡 点 ， 而 (1,2) 是 1- 型 
平衡 点 ， 其 部 分 不 稳定 流 形 收 敛 到 平衡 点 (0,0)。 对 系统 (11) 在 点 (1 2) 处 ， 由 于 


of 2— T2 =T} E 0 —1 
oz la -25 1l-m /ls, DU ec 0 
特征 值 Xi = V2 所 对 应 的 单位 特征 向 量 为 ai = (71/3, V/2/3)7, HIER Ao = -V2 所 对 应 的 单 
位 特征 向 量 为 as = (1/3,V2/3)7， 取 s = 0.1， 则 在 点 
(z, z) = (0.966667, 2.047133) 和  (z(9), (9) = (1.03333, 1.952867) 
处 用 波形 松弛 方法 对 系统 (11) 进行 求解 ， 其 雅 可 比 迭 代 格 式 为 
dg t» 


dt 


- 2o (70 H MUR 
(k+1) 
—-a(D-a( 9:9, ko012... 


画 出 轨 线 图 如 图 1 所 示 ， 其 中 虚线 为 用 波形 松弛 方法 计算 出 的 系统 (10) 渐 近 稳 定 平衡 
点 (0,0) 的 稳定 域 边界 流 形 ， 所 得 稳定 域 与 系统 实际 稳定 域 相同 。 
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图 1: 系统 (10) 渐 近 稳 定 平衡 点 (0,0) 的 稳定 域 


例 2 考虑 如 下 系统 88 


dzi 
ED 一 —T1 T I2; 
(12) 


d 
E = 0.12, — 2.025 — T? — 0.123, 
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Z1 一 2 
—— = T1 — 42, 
ds i 
T = -0.1zl + 2.022 + z2 + 0.123. 


系统 有 三 个 平衡 点 : (0,0) 是 系统 (12) 的 局 部 渐 近 稳定 平衡 点 ， (-2.55,-2.55) 是 系统 (12) 的 
1 型 平衡 点 ，(-7.45, -7.45) 是 系统 (12) 的 局 部 渐 近 稳定 平衡 点 。 对 于 系统 (13) 在 点 (-2.55,-2.55) 


9f 1 zi " 1 =i 
OT !(-2.55,-2.55) —0.14-2z,--0.322 2.0 -3.24925 2 ] 


特征 值 AM; = 3.370628 所 对 应 的 特征 向 量 为 ai = (1, 2.370628) ， 特 征 值 2 = —0.370628 所 对 
应 的 特征 向 量 为 aa = (1,0.629372)7， 取 < = 0.1， 则 在 局 
(x(x) = (—2.450000,—2.823063) 和 (zP, z$) = (—2.650000, —2.312937) 


(一 2.55, 一 2.55) 


处 ， 由 于 


处 用 波形 松弛 方法 对 系统 (13) 进行 求解 ， 其 雅 可 比 和 迭代 格式 为 


dz) — (een Q9 
dt i T2 >» 

deft” (k+1) (k+1) | (4 (0)? (k))3 

9» .17+tD) 2020? (337) +0.1(z1 ), k=0,1,2,... 


dt 
画 出 轨 线 图 如 图 2 所 示 ， 其 中 虚线 为 用 波形 松弛 方法 计算 出 的 系统 (12) 渐 近 稳 定 平衡 点 


的 (0,0) 和 (-7.45,-7.45) 的 稳定 域 边界 流 形 ， 所 得 稳定 域 与 系统 实际 稳定 域 相 同 。 
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图 2: 系统 (12) 渐 近 稳定 平衡 点 (0, 0) 和 (7.45, -7.45) 的 稳定 域 


例 3 考虑 Van der pol fT 7; £9! 
dzı 
— mi2; 
| (14) 


d 
audi = zı + A(z? = 1)z2, 


dzı = 
— = 12, 
| (15) 
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系统 仅 有 一 个 平衡 点 (0,0), 取 参 数 入 = 1, 在 (0,0) 较 远 处 取 一 点 (09, z0), 比如 (-1.8, 0.3)， 
从 此 点 开始 对 系统 (15) 用 波形 松弛 方法 求解 ， 其 雅 可 比 和 迭代 格式 为 


dz(*+1) _ 
dt a 

dig? (k) (k+l) k)A2, (k 

Tcu raf - aa, ko n2. 


画 出 轨 线 图 ， 可 得 系统 (14) 平衡 点 (0,0) 的 稳定 域 ， 并 且 可 得 系统 (14) 存在 不 稳定 极限 环 ， 其 
中 虚线 为 用 波形 松弛 方法 计算 出 的 系统 (14) 渐 近 稳定 平衡 点 (0,0) 的 稳定 域 边界 流 形 ， 所 得 稳 
定 域 与 系统 实际 稳定 域 相 同 ， 如 图 3 所 示 。 
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图 3: 入 = 二 1 时 系统 (14) 渐 近 稳 定 平衡 点 (0,0) 的 稳定 域 


5 小结 


我 们 用 波形 松弛 方法 计算 了 微分 动力 系统 渐 近 稳定 平衡 点 稳定 域 的 边界 流 形 ， 得 到 了 渐 近 
稳定 平衡 点 的 稳定 域 ， 由 于 已 经 得 到 微分 代数 系统 多 指标 85,16,18] 以 及 泛 函 微分 方程 nb9 波形 松 
弛 算法 的 收敛 性 结论 ， 我 们 可 以 尝试 用 波形 松弛 方法 计算 微分 代数 系统 ( 含 多 指标 ) 以 及 泛 函 微 
分 方程 的 稳定 域 或 者 收敛 域 问题 。 
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The Waveform Relaxation Method for Estimating the Asymptotic 
Stability Region of Nonlinear Differential Dynamic Systems 


LIN Xiao-lin 


(School of Science, Xi'an Jiaotong University, Xi'an 710049) 


Abstract: Estimating the asymptotic stability region of nonlinear differential dynamic systems is a 
problem applied in many fields. The paper outlines the computing methods for estimating the asymp- 
totic stability region of nonlinear differential dynamic systems. By analyzing the boundary manifolds 
of the asymptotic stability region, we present the idea of calculating the boundary manifolds of asymp- 
totic stability region with the waveform relaxation method and offer the algorithms. Firstly, find out 
all equilibrium points and determine which equilibrium point states in the boundary of asymptotic 
stability region of the stable equilibrium point. Secondly, use the negative system of nonlinear differ- 
ential dynamic systems to fix stable manifolds of equilibrium point setting in the boundary. Thirdly, 
the boundary of the asymptotic stability region of stable equilibrium point is composed of stable equi- 
librium points and their stable manifolds. Numerical examples further illustrate the correctness of the 
theoretical results in this paper. 

Keywords: nonlinear differential dynamic systems; region of asymptotic stability; waveform relaxation 
method 
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